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Abstract 

We study a model of densely packed self-avoiding loops on the annulus, related to the Temperley 
Lieb algebra with an extra idempotent boundary generator. Four different weights are given 
to the loops, depending on their homotopy class and whether they touch the outer rim of the 
annulus. When the weight of a contractible bulk loop x = q + q" 1 £ (—2,2], this model is 
conformally invariant for any real weight of the remaining three parameters. We classify the 
conformal boundary conditions and give exact expressions for the corresponding boundary scaling 
dimensions. The amplitudes with which the sectors with any prescribed number and types of 
non contractible loops appear in the full partition function Z are computed rigorously. Based on 
this, we write a number of identities involving Z which hold true for any finite size. When the 
weight of a contractible boundary loop y takes certain discrete values, y r = with r integer, 

other identities involving the standard characters K r , s of the Virasoro algebra are established. The 
connection with Dirichlet and Neumann boundary conditions in the 0(n) model is discussed in 
detail, and new scaling dimensions are derived. When q is a root of unity and y — y r , exact 
connections with the A m type RSOS model are made. These involve precise relations between 
the spectra of the loop and RSOS model transfer matrices, valid in finite size. Finally, the results 
where y — y r are related to the theory of Temperley Lieb cabling. 

SPhT-T06/155 

1 Introduction 



Boundary conformal field theories (CFT) have lately played an increasingly important role in statistical 
mechanics, condensed matter physics and string theory. In statistical mechanics, they appear in most 
probabilistic applications of geometrical models (see, e.g., [T] for a recent example), in particular 
through SLE 2]. In condensed matter, they contain all the information about fixed points in theories 
which are gapless in the bulk, such as Kondo systems or edge states in the fractional quantum Hall 
effect (see [3] for a review). In string theory, they provide for instance microscopic techniques to 
study D-branes in curved backgrounds [4]. On top of this, the study of boundary aspects is a crucial 
component of understanding and classifying CFTs at large [5] , and has lately played a crucial role in 
the solution of non rational CFTs [BJ. 

While progress in understanding conformal boundary conditions for rational CFTs has been consid- 
erable, the situation is not so satisfactory for non rational theories, which are however all too frequent 
in statistical mechanics applications. A case in point concerns the loop or cluster models, which — in 
one guise or another — are hidden behind most simple models of interest, such as Q-state Potts models, 
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0(n) models, RSOS models, polymers, and percolation. We are not aware of answers to most questions 
one might ask in this context, such as "what are all the conformal invariant boundary conditions", or 
"what are the partition functions for the O(n) models with Neumann boundary conditions", etc. The 
origin of this difficulty lies in our lack of understanding of the bulk CFTs, which exhibit non-rational, 
logarithmic features, and for which too little is known. Bulk exponents turned out to be tractable 
thanks to the Coulomb gas technique, but this technique has not been generalized to the boundary 
case with sufficient control yet (see [1] for recent progress in this direction). 

We put forward in this paper a proposal for what we believe are all the conformal boundary 
conditions of dense loop models. For each of those we determine the critical exponents, operator 
content and boundary partition functions, some of which have interesting probabilistic or geometrical 
interpretations. 

There seems to be much substance behind the results we uncover, and we hope to get back to 
the question in more details in the near future. In the present paper, we only present the leading 
arguments — which are based on previously published but unexploited results, as well as algebraic 
considerations — together with intensive numerical checks and some combinatorial proofs. 

To help the reader, we now give a quick summary of our results and notations. The boundary loop 
model (BLM) to be studied is defined on a tilted square lattice (see Fig. [3]), wrapped on an annulus of 
width N strands and circumference M lattice spacings. Loops cover all the edges, and interact in a 
specific way with the outer rim of the annulus, whereas they arc simply reflected by the inner rim (free 
boundary conditions). We denote by L the number of non contractible loops (note that L and N must 
have the same parity). Any loop has one of four weights (x, y, I or m, see Fig. |5}: I (resp. m) for a non 
contractible loop never touching (resp. touching at least once) the outer rim, and similarly x (resp. y) 
for contractible loops. We parametrize x = q + q^ 1 S (—2,2] by q = e l7T ^ p+1 " > (p real); the model is 
then critical with central charge (|2.7j) for any real values of y, I, m and is endowed with the U q (sl2) 
quantum group symmetry. We further parametrize y = y(r) as in (|2.4[) . Our central claim is that for 
any real r, and any L, there are two (distinct for L > 0) conformal boundary conditions: blobbed (resp. 
unblobbed) in which the outermost non contractible loop is required to (resp. required not to) touch the 
outer rim of the annulus. (When L = the two cases coincide.) The spectrum generating functions in 
these two cases are p.8[) . and the boundary conformal weights (critical exponents) h r , r ±L are read off 
from (|2.8|) . They combine to form the BLM partition function Z through the amplitudes (|3.1I|) . When 
p > I is integer, and when further r = 1,2, ... ,p the BLM model can be related to an RSOS model 
of the A p type with specific boundary conditions (three columns of fixed heights, see Fig. [T2")) through 
the rules (|4.4|) . In the latter case, Z can be written as a sum (j3. 18[) over irreducible representations of 
the Virasoro algebra. 

The paper is organized as follows. In Section [2] we review the algebraic framework used in our 
study (the blob algebra) along with a few key results. In Section [3] we define the BLM, classify its 
conformal boundary conditions, and give exact results for the associated critical exponents. Two 
appendices present a rigorous result on the amplitudes of the transfer matrix eigenvalues. This is 
used in Section 13.41 to write a number of exact identities — exact in finite size — relating Z to Virasoro 
characters K r s . In Section ^. 5l we discuss the case of Neumann boundary conditions for the loop model, 
identifying in particular the Neumann to Dirichlct boundary condition changing field. The relations 
to RSOS models are discussed in Section |U Finally, in Section [5) we comment on the relation between 
the blob algebra and the theory of Temperley Lieb cabling. Our conclusions — and the prospects for 
(much) further work — are given in Section [6l 



2 The blob algebra 

The Temperley Lieb (TL) algebra T/v (x) on N strands is defined by the generators (i = 1, 2, . . . , N—l) 
acting on strands i and i + 1 and satisfying the well-known relations 

e^ej = ejei for \i — j\ > 2 
e\ = xe, (2.1) 



2 



1 














X 




X 


< 


' 1 1 




1 




e 2 


e 3 


b 



Figure 1: Graphical representation of the action of the generators of the blob algebra, here shown for 
a system of N = 4 strands. 
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Figure 2: The word eieseee2eebeie3 = xyeie^e^ in the blob algebra B(x,y) on N = 7 strands. 

In [7] this was generalized into the two parameter "blob algebra" Bn(x,v), having an extra generator 
b, and satisfying in addition the relations 

b 2 = b 
eibei = yei 

etb = be, for i = 2, 3, . . . , N - 1 (2.2) 

For TJv(a;) it is well-known how to interpret these algebraic relations graphically in terms of the 
strands. The extra generator b marks the leftmost strand by adding a "blob" to it (see Figure [TJ. In 
this graphical representation any completed loop may be taken out and replaced by its corresponding 
weight (see Figure [5]). A loop with no blob gets the usual weight of x, while a loop with a blob gets a 
modified weight y. Note that several blobs on the same loop reduce to a single blob. Obviously, only 
loops touching the left border can be blobbed. 

This algebra has given rise to much work in recent years in the mathematical literature [jJJ. It has 
also been studied in the context of boundary conformal field theory [TT] with results that have some 
small overlap with ours. (In [TT] this algebra is called the "one boundary TL algebra" , a name we shall 
not adopt.) The blob algebra is in fact a quotient of the more general affine Hecke algebra (like TL 
itself is a quotient of the ordinary Hecke algebra). 

The representation theory of the two parameter algebra is richer than the representation theory 
of the TL algebra. For a given x = q + q , and assuming first that q = e* 7 is not a root of unity, 
exceptional cases occur whenever 

sin(r ± 1)7 . 
y = ; , r integer (2.3) 

sm r7 

[In the original paper [7] this corresponds to 77 = ^frj mod ir in the basic equation y = (q — q~ 1 e 2vq ^ / (l — e 2 ^) . 
In the case 

sin(r + 1)7 

y = ; , r integer (2-4) 

sinr7 

the spectrum of the hamiltonian 

H = — r — ( -ab - V e t j (2.5) 
tt sm 7 \ t^r J 

(the normalization guarantees unit sound velocity) has been studied in the continuum most recently 
in [TT] where it was found for any a > to give rise to the generating function of scaled gaps, in the 
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Figure 3: Vertices of the boundary loop model on the tilted square lattice. Bulk vertices are in any 
of two different states (corresponding to the generators 1 and e$). Left boundary vertices are always 
blobbed (generator b), and right boundary vertices are unblobbed (generator 1). 



sector with L non contractible lineal propagating: 



^Lo-c/24 _ g g 



Here, as usual, Lq is a Virasoro generator, 



c=l-—^— (2.7) 
P(P + 1) 



is the central charge for 7 = and 



[(p + ljr-^-l 

4p(p + l) V j 

are the conformal weights of the Kac table. Moreover, P(q) — Yi^=iO- ~Q n )- We stress that in (|2.6|) we 
have g = exp(27rir), where r is the standard modular parameter. As this meaning of q will be reserved 
for the argument of the spectrum generating functions, no confusion should arise with the other meaning 
of q as the quantum group deformation parameter q — e 11 appearing in the parameterization of x. 

The case r = 1 of (|2.6[) is the usual hamiltonian for the loop model with free boundary conditions, 
since in this case y — 2 cos 7 = x. In this case, the generating function of scaled gaps has indeed been 
known for a very long time [12] to be given by K\\+l. The case n > 1 as presented in [IT] appeared 
to be new, although it is in fact related to results in [12) (and has algebraic connotations in terms of 
representation theory of the blob versus the Virasoro algebra). In this paper, we shall discuss l|2.6p 
further, interpret it in the language of the loop model, and correct it whenever necessary. 

We stress that the independence upon a is a truly remarkable phenomenon: it can be interpreted 
by saying that once the algebraic structure of the hamiltonian is decided, the continuum limit does 
not depend on the (boundary) details. Another way to view this independence is the following. The 
R matrix of the loop model with spectral parameter u, acting on strands i and i + 1, can be written 
Ri(u) = 1 + f(u)ei, and f(u) = * s easily found by solving the Yang-Baxter equations. Writing 

similarly the boundary matrix as B(u) = 1 + g(u)b, the Sklyanin (or reflection) equations relate B(u) 
and Ri(u), giving rise to a solution for g(u) (see Eq. (40) in [T3]) that contains an arbitrary constant 
of separation (. The arbitrariness of ( is analogous to the a-independence discussed above. 



3 Boundary loop model 

We now want to study the blob algebra in the context of isotropic dense loop models, which are 
described by a transfer matrix instead of a hamiltonian. We recall that in the bulk, these models are 
defined by dense coverings of the (tilted) square lattice with self avoiding and mutually avoiding loops, 
each vertex allowing two possible configurations (see Figure [3] and also Figures fMHTSl below) . and each 
loop coming with a fugacity x. What happens at the boundary is the subject of this paper. 

1 A non contractible line is a strand propagating throughout the system. Thus, e; acting on two non contractible lines 
at i and i + 1 is zero by definition. Figure [2] has L = 1 non contractible line running from the bottom to the top of the 
figure (top and bottom are later wrapped onto an annulus). Note also that L and N must have the same parity. 
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Figure 4: (Color online). Numerical check of the conjecture h(y) — h r>r as a function of y, here for 
L = and p = 3 (the Ising model). The results were obtained from the finite- N corrections to the 
leading eigenvalue of the transfer matrix (|3.2p . obtained by exact diagonalization techniques, for widths 
up to N = 20. The left (resp. right) panel shows the choice A = 1 (resp. A = oo), but the extrapolated 
results (N — > oo) appear to be independent of A. The vertical lines represent the particular values 

map. 



The sum over TL generators in the hamiltonian is replaced in that case by a product 

L(iV-l)/2j [N/2] 

T = n (! + e 2i) II (i + ^-i) (3.1) 

i=l i=l 

where [■ ■ -J denotes the integer part. By analogy we will supplement this by boundary contributions 
so the full transfer matrix reads 

T = (1 + A6)T (3.2) 

In view of the a-independence of (|2.6[) . we would expect the critical exponents associated with 
T to be independent of A. This independence is checked numerically below (see Figure H]). From a 
geometrical point of view the most natural choice is then A = oo, so that after a trivial rescaling 
T = bTQ. This can be interpreted as a lattice model for which every loop touching the boundary gets 
a modified weight y instead of x. We will now study the conformal properties of this model (which we 
call the boundary loop model), as a function of x and y. 

3.1 Conformal boundary conditions and critical exponents 

The results from [12j and more recently suggest the following. Recall the parametrization x = 
2 cos 7; for each y(r) we solve (|2.4|) to get r. This gives a real number r £ (0,^)- The leading 
eigenvalue in the loop model transfer matrix with L non contractible lines should scale with conformal 
weight h(y) = h Ttr+ L- For when 7 = with p integer and r = 1,2, . . . ,p integer this is a rigorous 
consequence of jT^] and the loop-RSOS correspondence, as will be discussed in Section Q] below. 

Figure 0] serves the double purpose of checking this conjecture for L = and p = 3 (the Ising 
model) , and establishing the independence of the exponents on the choice of A in (|3.2[) . The agreement 
with the numerics is generally very good, and even in regions with strong corrections to scaling (in 
particular y = 0) it should be noted that the finite- N effects have consistently a trend and an amplitude 
compatible with the conjectured result in the thermodynamic limit. The choice A = 00 appears to 



5 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 2.2 2.4 



Figure 5: (Color online). Another check of the conjecture h(y) — h TtT for L = 0, here for p = 6 (the 
tricritical Potts model). 



minimize the amplitude of the finite-size effects, and accordingly we shall invariably adopt this choice 
for the subsequent numerical checks. 

Another check, still for L = but with a higher value p = 6 (the tricritical Potts model), is shown 
in Figure [5l 

From the point of view of the blob algebra, the parametrization (|2.4|) has nothing special compared 
with the other choice of sign in (|2.3|) . which we rewrite here as 

, sin(r / - 1)7 

w ) = ; ; (3-3) 

yv ' sinr'7 v ; 

Of course by using symmetries of the sine function we can write as well 

sin(p + 2 - r')-y 

V{r ] = sin(p+l-r07 (3 ' 4) 

so we see that the associated exponent reads as well h p +i- r > , p +i- r '+L = /v-i.r'-L, by the symmetry 
of (gH) . 

We now have to make things a little more precise and technical. The sector with L non contractiblc 
loops can in fact be considered from two points of view, depending on whether or not the leftmost 
non contractible loop is allowed to touch the left boundary. In Appendix A we discuss in details the 
sector structure of the transfer matrix VTq of (|3.1[) and show in particular that each of its eigenvalues 
corresponds to a definite choice: either the leftmost contractible loop is forced to touch the left boundary 
at least once, or it is forbidden from ever doing it. We shall henceforth refer to these two cases as the 
blobbed (resp. the unblobbed) sector. Note that in both sectors the contractible loops to the left of the 
leftmost non contractible loop may of course still touch the left boundary. 

We claim that these two cases both correspond to conformal boundary conditions, which are dif- 
ferent. For entropic reasons, the largest eigenvalue in the blobbed sector is obviously greater than the 
largest eigenvalue in the unblobbed sector. From the above arguments, the blobbed sector therefore 
has the exponent h r ^ r +L indeed. The unblobbed sector meanwhile has a different exponent /i r ,r-L- 
This can be checked numerically, and is illustrated in Figure |U Note of course that when L = the 
two results actually coincide, as they should, since the two sectors are then identical. 
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Figure 6: (Color online). Exponents h(y) in the presence of L = 2 non contractible loops, here for 
p = 3 (the Ising model). The lower (resp. upper) curves show the numerical results for the (constrained) 
boundary loop model. They agree with the conjecture h(y) = h r ^ r+ L (resp. h(y) — h r ^ r ^^) as claimed. 



3.2 Relation to the Potts model 

The boundary loop model is closely related to the Q — x 2 state Potts model at coupling (inverse 
temperature) J. Indeed, ignoring first boundary effects, the Potts model partition function can be 
written as [9] 

z= (e J -i) B Q c = Q v/2 Y,Q e/2 ( 3 - 5 ) 

clusters loops 

The first sum is over bond percolation clusters consisting of B bonds and C connected components. 
The second sum is over loops on the medial lattice that separate the clusters and their duals, with I 
being the number of loops and V the number of Potts spins. We have here supposed that the model 
is defined on a square lattice and stands at its critical temperature, e J — 1 = Q 1 ! 2 . The equivalence 
between the cluster and loop formulations is obtained by applying the Euler relation. Note that the 
local configurations of the loops correspond precisely to the first two vertices of Fig. [3] 

In the sector L = we have claimed that the exponent of the boundary loop model is h r , r in the 
parameterization (|2.4p . We now wish to check that this claim is consistent with known results on the 
Potts model. To that end, consider the Q-state Potts model in the same annular geometry as the 
boundary loop model. Denote by P± and P2 two points on the left boundary, and let boundary spins 
on the interval P\P 2 be constrained to take a subset of Q s states (with Q s < Q). In particular, when 
Q s = 1, this boundary condition corresponds to the Potts spins being fixed on the interval P\P 2 and 
free on the remainder of the boundary. The modified partition function reads 

Z( Pl P 2 )= £ (eS-l) B QC(^) = QV/ 2 Y Q i/2 (f) (3.6) 



clusters loops 



where we have used the Euler relation as before. The number of clusters (resp. loops) that touch P\P^ 
is denoted C{P 1 P 2 ) (resp. £(PiP 2 )), and obviously we have C(PiP 2 ) = t{PiP>z). We stress that C and 
I still denote the total number of clusters and loops. Now, (|3.6|) is a special case of the boundary loop 
model with the correspondence between weights 

Q = x 2 

Qs = xy (3.7) 
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In particular, for Q s = 1 we have r = p — 1 in (|2.4p . The corresponding special case of our general 
claim is therefore that the operator that changes the Potts model boundary conditions from free to 
fixed is ipp-i^p-i = <f>i,2- This indeed coincides with a well-known result of Cardy [10j . 

Another verification is furnished by Q = 3, Q s — 2. The claim is then that the operator that 
changes the boundary conditions from free to mixed is 02,2 with conformal weight ft.2,2 = 35 ■ This is 
again a well-known result. 



3.3 Spectrum generating functions 

Further numerical study shows that the spectrum generating functions for the blobbed and unblobbed 
sectors of the boundary loop model are simply the characters of generic irreducible representations of 
the Virasoro algebra, i.e., respectively, 

h r „r r — c/24 

Blobbed sector: Z% (r) = q ' t , , 

' (3.8) 



Unblobbed sector: Zi,{r) = — 



This can be related with the structure of the basis of the blob algebra and the absence of truncation of 
the Bratclli diagram [7] . The precise finite-size definition of Z*^ and Zl in terms of the transfer matrix 
blocks T£ and Tl (defined in Appendix A), is given in (|B.lj) - (|B.2[) : note that we have set j — L/2 in 
Appendix B. 

The leading behavior in these expressions defines the exponents h r<r ±L and has already been checked 
in Figs. above. The coefficients of the terms up to level 6 in the development 

1/P(q) = 1 + q + 2q 2 + 3q 3 + 5q 4 + 7q 5 + llq 6 + ... (3.9) 

have been verified by computing the first 32 eigenvalues of T, for sizes up to N — 24, and looking for 
integer gaps in the spectrum of critical exponents. For definiteness we have concentrated on the case 
p = 3 and L = 2. Independent computations were made in the blobbed and the unblobbed sectors. 
Moreover, the verification was made both for y — x and for a generic value of y, and in cither case the 
absence of singular vectors up to and including level 6 was ascertained. 

It is important to stress that the spectrum generating functions are not given by (|2.6p for the 
BLM. The full loop transfer matrix actually contains more information, but can be truncated when r 
is integer. 



3.4 Partition function identities 

We are now interested in the situation where the non contractible loops wrap around the annulus. The 
question then arises of which weight should be given to these loops (while our results for the exponents 
have some overlap with [llj . the following has never appeared before). We will give in general the 
weight I to non contractible loops that do not touch the boundary (i.e., the outer rim of the annulus), 
and weight m to those that do (there is at most one). We now claim that the full partition functions 
are given by adding sectors Z^(r) and Zl(t) with the following amplitudes. We parametrize I, m in 
terms of two numbers a, (3: 

I 

m 

Then the amplitudes are 

D* L 
D L 



2 cosh a 
sinh(a + (3) 



sinh/3 




smh(La + 


P) 


sinh (3 




sinh(La — 





sinh(-/3) 



(3.10) 



(3.11) 
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Figure 7: Pictorial representation of the recursion relation (|3.13|) . 



A complete proof of this statement is presented in Appendix B. Alternatively, we can argue that, 
because one can get from a sector L to a sector L + 2 by adding two non contractible loops at the right 
of the annulus (here seen as a periodic strip), which can or not get contracted with the first L ones, 
the generic amplitudes have to obey a recursion relation (which has a deep algebraic nature, see [H]). 
Considering first the amplitudes D° L for the simpler problem with transfer matrix To (i.e., without the 
boundary generator b) one has 

D° L D° 2 = D° L + D° L+2 + D° L _ 2 (3.12) 

with initial values Dq = 1 and D 2 = I 2 — 1 . A pictorial rendering of this relation is shown in Figure [7] 
The solution reads explicitly D° L = Ul(1/2), where Ul is the Tth Chebyshev polynomial of the second 
kind. Turning now to the boundary loop model we must have similarly 



D L D° 2 = D L 



D* L D° 2 



I), 



D L+ 2 
D*L+2 



L-2 



L-2 



(3.13) 



The initial values can then be determined to be Dq = Dq = 1, D 2 = I 2 — Im — 1 and D 2 = Im — 1, 
from which the general formulas (|3.11[) follow. In addition to the actual proof of Appendix B, we have 
checked (13. 1 1[) by formal manipulations of transfer matrices up to size N = 6, along the lines of [15] . 

The amplitudes can then be used to write down the general partition functions in the case N even 
(we set L = 2j) 



Z 



-c/24 



sinh(2ja + ,9) q h 

3=0 



sinh (3 



P(q) 



sinh(2jo! — 0) q 



sinh (3 



(3.14) 



This partition function correspond to the most general case represented in Figure [5] contractible loops 
in the bulk get the weight x, those touching the boundary a weight y, non contractible loops not 
touching the boundary a weight I and the others a weight m. 
The simplest situation occurs when 



rn 



sinh(7' + l)a 
sinh ra 



(3.15) 



that is, the parameter m has the same formal algebraic relationship with / that y does with x. Then 
(3 = ra — the same r as in (|2.4[) — and one can write 



Z 



-c/24 



E 

3=0 



sinh(2j + r)a q hr 



y 

sinhra P{q) ~[ sinhra P{q) 
If moreover one is in a degenerate case where r is integer, one can reorganize the sum (|3.16[) into 
sinh(2j + r)a 



sinh(2j — r)a q hr ^- 2 i 



(3.16) 



3=0 



sinh i 



K r 



~+2j — q 



-c/24 sinh ( 2 J - r ) a q K - r - 2 ' 



3=1 



sinh / 



P(l) 



(3.17) 
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Figure 8: (Color online). The most general model studied in this paper. We distinguish four types of 
self-avoiding loops on the annulus: contractible loops touching (blue color in the figure) or not touching 
(green) the outer rim of the annulus, and similarly for non contractible loops (red or purple). 



By pairing up terms with j and r — j in the second sum, this can in turn be rewritten as 

Z= ± ^h(2i + r). 

smhra 

j=-[r/2] 

where [• ■ -J denotes the integer part. When r is even, the contribution from j = —r/2 actually 
disappears. We thus get a sum over irreducible representations of the Virasoro algebra. 

We claim that the subtractions occurring in the partition function do occur in finite size as well. 
These subtractions involve the conformal weights h rjr+ 2j and /i r ,-r-2j and correspond respectively to 
the blobbed sector with L = 2j non contractible loops and the unblobbed sector with L = 2j + 2r. In 
other words, we claim there are level coincidences in finite size between the blobbed and unblobbed 
sectors when r is an integer: this in fact follows from the theory of representations of the blob algebra. 

To be more precise, one can make sense of expressions such as (j3.18[) in finite size by replacing the 
definition (|2.6p of the characters K r , s by traces of transfer matrix blocks, as outlined in Appendix B. 
Care should then be taken that the annulus is wide enough to accommodate the prescribed number of 
non contractible lines, which amounts to replacing the upper limit in the summation (13.18|) by N '/2. A 
numerical check of the level coincidences is shown in Tablc[TJ for the case N — 6, p = 6 and r = 1 (i.e., 
y = x). We find indeed that the level spectrum of the transfer matrix (see Appendix A) block Tij+ir 
is a proper subset of that of T 2 *j, f° r all allowed values of j (i.e., j — 0,1, ... , N/2 — 2r — 1). We have 
checked this statement for several other values of N, p and r, including on examples involving many 
more levels. 

Although we have restricted to ./V even, similar results can be written for N odd, with the difference 
that L is odd (and thus never vanishes). The generating function then reads 



q -c/2A 



^ S inh[(2j + l)a + j3] q h -.-+^+ 1 ^ s inh[(2j + I) a - 0\ q h ^-sj- 



sinh/3 P(q) ^ sinh/3 P(q) 



(3.19) 



3.5 Neumann boundary conditions 

If we go back to the derivation of the partition function for the 0(n) model [16] we see that Dirichlet 
boundary conditions for the 0{n) variable will translate into having a loop extremity on every point 
of the boundary, and the open lines thus obtained get a weight one (since the 0(n) variable is fixed to 
say S = (1, 0, ... 0) on the boundary.) Let us assume this carries over to the fully packed case and our 
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Table 1: Complete set of levels for the various transfer matrix blocks (cf. Appendix A) in the case of 
N = 6 strands, for p = 6 (the tricritical Potts model) and r = 1 (i.e., y = x). The left column shows 
— A -1 log A (rounded to 12 digits), where A > 1 is the transfer matrix eigenvalue, and the remaining 
columns show the multiplicities of A within the various blocks. A blank entry denotes zero multiplicity. 
The level coincidences mentioned in the text are clearly observed. (There is one extra coincidence of 
the level -0.085859268861, between T 2 * and T|.) 
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Figure 9: (Color online). Configuration in which the two boundaries of the annulus (here shown as a 
periodic strip) have respectively Dirichlet and Neumann boundary conditions. RSOS faces are shown 
in black, and the loops that separate clusters at constant height are green. The open lines going from 
the left to the left boundary carry unit weight. 



geometry (see Figure [9]) , where we thus demand that every point on the boundary looks like the top 
diagram on Figure 1 1 01 which we call a fork. We also require that open loops thus formed all carry a 
weight unity. It is then easy to see that this is equivalent to marking loops touching the boundary with 
a blob having parameter y = 1. In turn, with the usual parametrization for x we get the associated 
conformal weight to be 

/iTwist = h p/2tP/2 = 16p (p +1 ) (3- 2 °) 

since y = 1 in (12. 4p when r = p/2. This corresponds to a twist operator, or the dimension of the 
boundary field changing boundary conditions from Neumann to Dirichlet. 

We note that this dimension of the twist operator agrees with the one proposed in [17] after some 
reinterpretation of the results. The formula given in this reference is slightly different 

krwtat-fcaJMJl- 16p(p+1) ( 3 - 21 ) 

but turns out to hold for the dilute phase of the 0{n) model [TB] only [TB]. For the dense case, it has 
to be replaced by (|3.20|) . 

We can now write the partition function with Dirichlet boundary conditions on one rim of the 
annulus and Neumann on the other, simply by setting r — | in our formulas. An interesting limit to 
consider is then p — > oo, m — 1, 1 = 2 where we find from (|3.14|) 

Z = £ <L_^ ( 3. 22 ) 

j=— oo 

which is the usual Dirichlct-Ncumann partition function for the free boson. The presence of the two 
kinds of terms h r ,r±L in (|3.14j) is crucial to recover this limit. 

Note that in the limit p — > oo, the exponent h(y) becomes a step function: 

(J ^T~ for y < 1 
hi !n ={ Q^t fory=l (3.23) 

TT for y > 1 
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Figure 10: (Color online.) Two possible interactions between the loops and the left boundary: fork 
(top image) and blob (bottom image). 
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Figure 11: (Color online). Exponents h(y) with L = non contractible loops, in the limit p — > oo. The 
exact result is a step function, passing through the value h(l) = (shown as a big cross). 
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Figure 12: RSOS model in an annular geometry, for N = 6. Heights on the rightmost and the two 
leftmost layers are fixed as shown. Time flows vertically, and there are periodic boundary conditions 
in the time direction. 

The numerical check, shown in Figure 1111 is compatible with this behavior, although the finite-size 
corrections are of course large near the step. 

Note that in general blob and fork do not coincide: this is true only when the associated weights 
are equal to unity. One can play the game of defining a fork algebra such that the top diagram in 
Figure ITD1 defines the fork operator /. Required relations to give to every open loop a weight z are 
then obviously 

f = zf 

d/ei = zei (3.24) 
But, by a rescaling / — > f/z we get the same relations as the blob relations with y = 1! 

4 Relation to RSOS models 

We now want to tackle the boundary loop model by another route. Recall that in the numerical studies 
of Saleur and Bauer [T2] it was found that for A p RSOS models [with central charge (|2.7p ] the annulus 
partition function is exactly the character \da when the following boundary conditions are imposed: 
all heights on the right boundary of the annulus are fixed to a (Dirichlet boundary conditions) while 
on the left boundary, the heights on the boundary are fixed to 6 and those in the layer next to the 
boundary are fixed to c, with d = min(6, c)@ Note that without the constraint on the next-to-leftmost 
heights, each height in that layer may take either of the values 6+1 and 6—1. This is illustrated in 
Figure Q2J 

Let us now see how this choice of boundary conditions translates in the loop model. To do so, we 
first note that if the loop model is defined on N strands, a time slice of the RSOS model is defined by 
N + 1 heights. A state of the model is a collection of those heights, denoted \h,h, ■ ■ ■ , In+i), and we 
have the RSOS constraint Zj = 1,2, ... ,p with \h+i — k\ — 1. Figure [T2l corresponds to a case of even 
N. We will restrict to this case to start. We will also restrict to the case where c = &+ lso6 = tf, 
c = d + 1; note that then d — a and N have the same parity. 

We next recall briefly how the Temperley Lieb generators act in the RSOS representation. Taking 

2 In order to respect the use of a, b and c in Ref. [121 . note that throughout Section [4] a is not the parameter of 112.51 1. 
b is not the blob generator in H2.2I I. and c is not the central charge 112.71 1. 
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Figure 13: Labelling of RSOS heights around a lattice face. 



the time to flow downwards in Figure 1131 the generator acts as 

{SiA') 1/2 

Bi . . . , k-i,li, li+l, . . . , ^W+l) — 5(k-i, U+l) 2, FT^ I' 




. . . , li+i, ■ ■ ■ , In+i) (4.1) 



where the Si are the components of the Perron-Frobenius eigenvector of the adjacency matrix of A 
and read explicitly 



where we have introduced the ^-deformed numbers [l] q = (q — q~ l )/(q — q 1 )- 

The transfer matrix of the RSOS model has the usual form (|3.ip . but now in terms of the e$ defined 
by (|4.ip . The graphical expansion of the partition function is obtained by taking, for each elementary 
face transfer matrix 1 + e*, either the identity — in which case a vertical bar is drawn diagonally across 
the face, indicating that It = l[ for this term — or the Temperley Lieb generator ej, in which case a 
horizontal bar is drawn, indicating that = U+i. One gets in this way clusters of constant heights, 
which can be separated by non intersecting loops drawn on the dual lattice. 

The analysis of the weights needs some modification as compared to Pasquier's original treatment 
[T9] of the RSOS model on a torus, but some of the key elements can be taken over. We first consider 
the case of li = I2N+1 = d; the heights I2 are also fixed and we shall assume I2 = d + 1 (the case 
I2 = d — 1 being similar). We further assume that there is at least one cluster connecting the left and 
right boundary of the annulus (in the limit where the aspect ratio p — M/N — > 00 this is almost surely 
the case). There is then L = non contractible lines. Since all clusters connecting the two boundaries 
have the same height d, we might just as well identify them. In the graphical expansion this can be 
represented by adding extra vertical bars on the two boundaries. 

With this identification it follows that any loop is at the junction between exactly two distinct 
clusters. Now orient every loop in the clockwise direction. When traversing any loop along this 
direction, the cluster to its right is said to be surrounded by the loop, whereas the cluster to its left 
is said to surround the loop. Let us now represent each distinct cluster by a node, and each loop by 
a directed link going from the cluster that it surrounds towards the cluster that it is surrounded by. 
This then defines a directed rooted tree, where every link is oriented towards the root, the root being 
the unique cluster connecting the two boundaries. 

The weight appearing in (|4.1[) is then distributed on individual loop turns as follows: Consider a 
loop that surrounds a cluster at height I and that is surrounded by a cluster at height k. When making 
a right (resp. left) turn whilst being tangential to a horizontal bar, the loop picks up a factor (Si/Sk) 1 ^ 2 
(resp. (Si/Sk)^ 1 / 2 )- Turns being tangential to a vertical bar do not carry any weight. The complete 
loop, being clockwise, makes four more right than left turns, but only two of those excess turns carry 
any weight, so the complete weight is Si/Sk- 

The directed tree is now undone by summing over the heights of its nodes. We start at a leave node. 
For any fixed height k of the node adjacent to the leaf, the height of the leaf node can be I = k ± 1. 
Summing over / gives a weight|f| (Sk+i + Sk-i)/Sk — 2 cos f^rj — note that this is independent of 

3 This is true even when I = 1 or I = p, since Sq = Sp+i = from (14.21 ). 




(4.2) 
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Figure 14: (Color online.) Configuration of the RSOS clusters (in red color) for a case where a = d. 
Periodic boundary conditions are imposed in the time (vertical) direction. 



k — which can be attributed to the loop represented by the link directed away from the leaf. We then 
remove the leaf and its outgoing link, and proceed iteratively, moving always from the leaves and 
towards the root. A special case occurs when summing over the heights of those nodes adjacent to the 
root whose outgoing link corresponds to a loop that touches the left boundary. Indeed, the heights of 
those nodes have been fixed to d + 1, and we obtain then the weight Sd+i/Sd = [d+ l] q /[d] q . Finally, 
when the whole tree has been undone, only the root node remains, but since its height is fixed it 
contributes no additional weight. 

Consider now the loop model with weights x = 2 cos f^frj arL d V = [d+ l] q /[d] q — i.e., setting r = d 
in (|2.4p — and no non contractible lines, L = 0. Its configurations are in one-to-one correspondence 
with those of the RSOS model treated above and the weights x, y are the same. The spectrum of Ti oop 
therefore contains the complete spectrum of TrsoSj and we have verified this numerically. In particular 
the leading eigenvalues of these transfer matrices coincide. However, Xi oop also contains eigenvalues 
not present in Trsos- This was to be expected, since the loop model contains non-local information 
not present in the RSOS model (allowing in particular its definition for non-integer p). 

An example of the loop-RSOS equivalence is illustrated in Figure [T4j still for the case d = a. 
Denoting by l\, . . . , k the heights of the five cluster^ (labelled from top to bottom) we find the overall 
weight (for simplicity we denote [l] q = I): 

(kk) 1/2 (hh) 1/2 (hk) 1/2 (kk) 1/2 {hh) 1/2 k = k k (43) 
h k k k k k k k 

Note that k = k by the periodic boundary conditions in the time direction. Each of the two factors 
on the right-hand side is associated with a loop. Since k can take values d ± 1 while k = k — d and 
?4 = li = d + 1, the overall weight is 2 cos ( ^pj^ x = xy as claimed. 

While this construction clearly works when there are no non contractible lines running through 
the annulus (ie, a = d), things are more delicate when there are. Figure [T51 shows for instance that a 
configuration with two non contractible lines does not necessarily contribute to d = 3, a = 1. A little 
thought suggests to simply eliminate non contractible loops touching the boundary in this case. 

We have indeed checked numerically that when d — a is negative, the eigenvalues of the RSOS model 
are all found in the loop transfer matrix for the blobbed sector, and the leading eigenvalues coincide. 
Meanwhile if d — a is positive, this is true provided one considers instead the unblobbed sector for the 
loop model. In both cases one needs to have L = d — a and y = y(d). 

4 Note that in this example the subscripts arc used differently than in 1 14. It . 
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Figure 15: (Color online.) The choice a = 3, b = 4, c = 1 is not compatible with the configuration 
of loops shown on the left panel (even though there are two non contractible loops). It is however 
compatible with the configuration shown on the right panel. 

This corresponds simply to accepting configurations such as the one on the right of Figure (|15[) . 
but not the one on the left. 

Of course the RSOS model contains more situations: one can decide to have c = b — 1 instead, or 
to have N odd. The complete set of rules is as follows: 

d — a > and has same parity as N: L = \d — 

d — a < and has same parity as N: L = \d — 

d — a > —1 and has opposite parity of N: L = \d + 

d — a < — 1 and has opposite parity of N: L = \d + 

The first two cases coincide with what we discussed if iV is even and generalize it easily if N is odd. 
In these two cases, the left boundary necessarily sees b = d + 1, c = d, and the value of y follows from 
the analysis of the loop model. 

The last two cases meanwhile require that b = d, c = d + 1. The new value of y giving the correct 
weights to the boundary loops is thus y = [d] q /[d+l] q , which coincides with y = [p + 1 — d] q /\p — d] q . 

Note that the rules (|4.4|) are compatible with the global symmetry Zj — > p — k of the RSOS model 
configurations and weights (|4.2[) . On the level of the sectors this reads {d,a) — > (p — d,p + 1, a) and is 
nothing but the usual symmetry of the Kac table (|2.8[) . 

An extensive numerical check of (|4.4[) is given in Table [51 We show all levels observed for the 
RSOS model with N = 6 and p — 6, along with the sector label (d,a), and we give the corresponding 
multiplicities of these levels in the loop model, where each sector Tj*'(r) is characterized by the number 
of non contractible lines L, the blobbing of the leftmost string (* denotes the blobbed sector), as well 
as r which determines the boundary weight y(r) through (|2.4|) . Remarkably, all of the RSOS levels 
are also observed in the loop model, with the sector given precisely by the rules (|4.4p . This observation 
extends to other values of N and p (with any parities) . 

Note that the dimension of the loop (resp. RSOS) model transfer matrix does not (resp. does) 
depend on r. This apparent paradox is resolved by the fact that a given loop model sector contains in 
general extra eigenvalues (not shown in Table [2]) which are not present in the RSOS model. It should 
also be noted that the dominant eigenvalue in a given RSOS sector is always observed to be dominant 
as well in the corresponding loop model sector (|4.4|) . 

Finally, we should mention a couple of fine details. In general the levels are not observed to be 
degenerate. One exception visible in Table [2 is that the level / = —0.036787385047 is present in two 
distinct sectors, but since this is so both on the RSOS and the loop side, definite labels respecting (|4.4[) 
can be assigned as shown in the table. Another exception is the level / = —0.085859268861 which 



a\, y = y(d), unblobbed 

a\, y = y(d), blobbed 

1 — a\, y = y(p — d), blobbed ' 

1 — a\, y = y(p — d), unblobbed 
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Table 2: Complete set of levels for the various sectors (d,a) of the RSOS model, here for size N = 6 
and p = 6 (the tricritical Potts model), cf. Table[TJ The left column shows —TV -1 log A (rounded to 12 
digits), where A > 1 is the transfer matrix eigenvalue. The remaining columns show the corresponding 
multiplicities of A within the various sectors T£ of the loop model (the asterisk * denotes the blobbed 
sector, see Appendix A). A blank entry denotes zero multiplicity. The boundary weight y = y(r) is 
given by (|2.4p . with r indicated in the top of the table. 
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apart from its assignment to T4 (r = 1) by (|4.4[) appears also in two more loop sectors, T£(r = 1) and 
T^(r = 1). This type of degeneracies can be explained by quantum group arguments. 

Recall now that the character of the irreducible representation of the Virasoro algebra with highest 
weight hda can be written as [20] 

00 00 

Xda = K d,a+2(p+l)n ~ ^ K cL2n (p+l)-a (4-5) 

n— n— 1 

We wish to recover this character from the knowledge of the loop model partition function. This involves 
as usual an infinite series of additions and subtractions of sectors which is made possible in finite size 
by the quantum group symmetry. We suppose we are in the situation with heights b = d, c = d + 1. 
To go from d on the left side to a on the right side we need to sandwich between the left and the right 
the adjacency matrix of the A p diagrams so that a random walk on this diagram, from boundary to 
boundary, hopping from non contractible cluster to non contractible cluster, takes one from b to a. All 
the steps are identical with those in [12j where partition functions with boundary conditions fixed at 
the leftmost and rightmost sides were computed. Using the eigenvectors of the adjacency matrix as in 
equation (4.15-4.19) of [12], the required expression is thus (we have set H = p + 1) 



2 ^ . / irpd\ . firpa\ ( inp N 

h ^ sm hrj sm hrJ z [ a = if 13 = da . 



(4.6) 



Here Z is calculated with a running over and for each a the value of m — the weight of non 
contractible loops — follows from the same argument as in the case of the leading exponent r = 1 in the 
bulk case, implying P — da so m — s ^^° , a — Using the expression (|3 . 1 8[) we see that the 
functions Kd y d+2j (r = d) get a combinatorial factor 



Jj \ [ cos ^' + d-a)^- cos(2j + d + a)^\ (4.7) 

Since d — a and 2j have the same parity, this select the conditions 

2j + d - a = 2niH 

2j + d + a = 2n 2 H (4.8) 

for the first and second term respectively. Remembering that j runs only from — [r/2] to infinity, 
and that when the foregoing conditions are satisfied one gets a factor from the sum cancelling the 
prefactor in (|4.6p . it follows that 

oo oo 

Z = ^ K d>a+ 2nH — K<i,-a+2nH (4-9) 
n— n— 1 

which is the result we wanted. 



5 More algebraic considerations 

There are many ways to think of the blob algebra. We would like now to think of it within the 
theory of cabling, i.e., tensor products of spin 1/2 representations of U q {sl2)- Consider therefore 
r + 1 representations of spin 1/2, and suppose we wish to project them on the maximally g-symmetric 
representation. The object doing this is the g-symmetrizer, whose expression is well known by induction 
to be [21] 

(S r _(-i(ei, . . . , e r ) — S r — t r S r e r S r (5-1) 
with the boundary condition S\ = 1. Here the t r are numbers given by 

tr = . f 11 ^, (5.2) 
sm(r + 1)7 
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Figure 16: (Color online.) Ghost strings 2 — r, ... ,0, represented by thin red lines. They are sym- 
metrized among themselves and with the first real line (the symmetrization is shown symbolically by 
horizontal green bars). 



and the are TL generators which obey (12. 1|) as before and act on the tensor product of the i'th and 
(i + l)'th spin 1/2 representation. These generators will not be identified with the ones used in the 
previous sections however (see below). The first values are well known: S± = 1, S 2 = 1 — 2 cos ' e ^ c - 

Let us now add to our system r — 1 "ghost strings" as in figure [TBI (labelled 2 — r, . . . , 0) on which 
generators e2- r , ■ ■ ■ > eo act, and let us symmetrize on these ghost strings and the first one of our system. 
We then define b r through 

b r = S r (e 2 - r , ■ ■ ■ , e ) (5.3) 

We have then 

h = 1 

b 2 = 1-^- (5.4) 
2 cos 7 

63 = 1 -t 2 (e_i +e ) +tit 2 (e-ie + e e_i) 
and so on. By construction, the 6 r 's are projectors 

b 2 r = b r (5.5) 

It is easy to show that they satisfy moreover 

sin(r + l) 7 

e\b r e\ — : e\Jr r —\ (,5.o) 

smr7 

where we denote by P r -i the symmetrizer on the ghost strings only, P r _i = iS r _i(e2- r) ■ ■ • , e_i). Of 
course, [ei, P r -i] = for all i > 1. We can thus consider a modified version of the TL algebra where 
instead of the generators we have the modified generators Sj = P r -\ei. They obviously satisfy the 
bulk TL relations (|2.1jl . since P r _i is a projector. Moreover we now have the relations l|2.2[l of the 
blob algebra, with b = b r and (|2.4jl . We have thus made the link with the foregoing discussion and 
shown that for r an integer, the boundary conditions corresponding to the value (|2.4[) can be obtained 
by adding ghost strings on the left boundary and symmetrizing them with the first "real string" in the 
system. 

This should not come as a surprise: it is easy to show using Pasquier's 6j calculations [H] that the 
insertion of b r through the ghost strings construction translates into RSOS language by having heights 
increasing linearly from the left hand side of the ghost strings up to the first height from I = 1 to 
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I = r, and then because b r symmetrizes on r strings and thus on the first real string as well, from I = r 
to I = r + 1. This proves that the transfer matrix bTo describes exactly the mixed RSOS boundary 
conditions. 

Note that one does not have to introduce the ghost strings. Algebraically, the same would be ob- 
tained by replacing these r— 1 strings by a spin (in which case, P r -i = 1 identically). The operator 
b r then amounts to projecting the product of this representation and the first spin 1/2 representation 
onto spin r/2. 

In this form, the identification was already mentioned in |12j . 

The formulas obtained in this section should match the ones which appeared in a recent paper by 
Pearce et al. [23]; the derivation there is based on boundary integrability, and does not refer to the 
blob or boundary Temperley Lieb algebra, as far as we can see. 

6 Conclusion 

In conclusion, it is important to stress that we have found a continuum of conformal boundary con- 
ditions for the dense loop models. How to incorporate them into a consistent conformal field theory 
remains an open problem. Note that none of these boundary conditions involves the number of times 
loops touch the boundary, which would correspond somehow to modified spin-spin couplings on the 
boundary. Rather, giving a different weight y to loops touching the boundary can be interpreted in 
the 0(n) model most easily as restricting the degrees of freedom on the boundary to take values in a 
sub manifold of 'dimension' y. This is not without reminding us of results in the WZW cases [2~i] . 

The results described in this paper point to many further directions. Among those are: 

- geometric applications of the generating functions on the annulus 

- derivation of the Bethe ansatz equations and of the spectrum of scaled gaps 

- extensions to the dilute case 

- extension to the double boundary case 

- study of boundary conditions in c = 1 theories 

We hope to report on these soon. To conclude, we now give one example of application. We consider 
the case x — 0, c = —2 and the partition function with x = I = 0, y = m = 1 (i.e,. loops touching 
the boundary get a weight one, whether contractible or not, the others are not allowed). After some 
simple manipulations, ()3 . 14(1 can be written as 

— 1/24 00 

z = Vt £ (-i)V 4 >- 1)2/32 (6.1) 

j=-oo 

This is easy to interpret geometrically as the partition function of a gas of dense loops, all contractible, 
which are all constrained to touch the left boundary. The loops can also be replaced by trees. Mean- 
while, this partition function can also be interpreted in the symplectic fermion theory |25j . 

Acknowledgments: H. Saleur thanks I. Kostov and V. Schomerus for many interesting discussions, 
and for communication of their unpublished results [181 125j . 

Note added in proof: After the completion of this work, I. Kostov has studied (in the preprint 
|hep-t h/070322i"| the coupling of our boundary loop model to two-dimensional quantum gravity. His 
results corroborate those presented here. 

A Transfer matrix structure 

In this Appendix we discuss the construction and structure of the transfer matrix of the boundary loop 
model, corresponding to taking the limit A — > oo in (|3.2p . 
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Figure 17: Ordering of the states of the transfer matrix T, here for width N = 4 strands. 



We recall the most general case of the boundary loop model, cf. Fig. [H The model is defined on 
an annulus with the exterior boundary being distinguished. Each loop touching at least once this 
boundary gets a weight y if it is contractible (i.e., homotopic to a point), and m if it is not. Each loop 
that never touches the distinguished boundary gets a weight x if it is contractible, and I if it is not. 
For an annulus of width N strands (we assume N even) and circumference M, this defines a partition 
function Zpf y M{x,y,l,m) which can be expressed in terms of the M'th power of the transfer matrix 
T = bTo, where T is given in terms of TL generators by (|3.ip . and the generators b, ej obey the blob 
algebra B N (x,y) defined by (|2~Tj) ~ (|2~2"l) . 

For simplicity we henceforth represent the annulus as a rectangle of width N and height M , with 
periodic boundary conditions identifying its top and bottom sides. The distinguished boundary is 
taken to be the left side. The transfer matrix then acts on states which can be depicted graphically 
as fully-packed non-crossing link patterns within a slab bordered by two horizontal rows, each of N 
points. A link joining the top and the bottom row is called a string, and any other link is called an 
arc. The action of a word in Bi<i{x,y) on a state is obtained by adjoining the word to the top row of 
the slab (i.e., time propagates upwards). The 20 possible states for N = 4 are represented in Fig. [T7] 

Note that links touching a point on the left boundary are necessarily blobbed (shown by a circle in 
Fig. [TTf . Only links up to and including the leftmost string can be blobbed. The states can be ordered 
as follows: First we sort the states according to a decreasing number of strings L. For fixed L > 0, we 
place first the states in which the leftmost string is unblobbed. And finally we group together states 
with fixed L and fixed blobbing of the leftmost string, according to the link/arc configuration of the 
lower row of the slab. This gives the order of the rows of Fig. [T71 The ordering of states within each 
row is according to the link/arc configuration of the upper row of the slab. 

With this ordering of the states, T has a blockwise lower triagonal structure, with each block 
corresponding to a group of states as defined above. The reason is that acting by e; can annihilate 
two strings (if their positions on the top side of the slab are i and i + 1) but cannot create any strings. 
Likewise, acting by b can blob the leftmost string, but it cannot subsequently be unblobbed. The 
triagonal structure implies that the eigenvalues of T are the union of eigenvalues of the blocks on its 
diagonal. Moreover, blocks differing only by the configuration of the bottom row are identical. For 
the purpose of studying only the spectrum of T the bottom row can therefore be completely forgotten, 
leading to a much smaller transfer matrix. 
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Figure 18: Reduced states for N = 4. 

Indeed let us define a reduced state as a non-crossing link pattern on N points. A (full) state can 
be turned into a pair of reduced states by cutting each of its strings and pulling apart the upper and 
lower parts. For convenience, a cut string will still be called a string with respect to the reduced state. 
For N = 4 there are 8 reduced states, shown in Fig. [TSl 

The blocks on the diagonal of T are denoted Tl and T£ , where L is the number of strings, and the 
presence (resp. absence) of the asterisk (*) indicates that the leftmost link is blobbed (resp. unblobbed). 
Note that the blocks Tl and T£ can be constructed in terms of the reduced states. The numerical 
studies of the spectral properties of T reported in this paper were done by diagonalizing these blocks 
in the basis of reduced states. Their dimensions read (see Appendix B for proofs of closely related 
statements) 

dimTi = (, N ~ 1 ,, | fori = 0,2,...,JV-2 
\{N - L-2)/2j ' ' ' 

dimT£ = ( ",,,_) for £ = 2,4, ...,N (A.l) 



\(N-L)/2 

With the terminology being fixed, the annulus partition function can now be written as [15] 

Z N , M (x,y,l,m) = (u\T M \v) . (A.2) 

Here, the right vector \v) is the unit vector corresponding to the state with N strings. The left vector 
(u\ is obtained by identifying the top and bottom rows for each state; counting the number of loops of 
each type (contractible or not, blobbed or not) gives the corresponding weight as a monomial in x, y, 
I, and m. For instance, with N = 4 and the ordering of the states shown in Fig. 1171 we have 

\u) = (1,0,0,. ..,0) 

\v) = (l 3 m,l 2 y,lmx,lm,y, . . . ,xy) (A. 3) 

B Exact eigenvalue amplitudes 

The goal of this Appendix is to provide a rigorous combinatorial proof of the amplitude formulae (|3.11[) 
by generalizing the working of [26] . The discussion assumes knowledge of the transfer matrix blocks 
Tl and T£ defined in Appendix A. 

Following [26] , we introduce the characters 

K k = Tr {T 2k ) M , K* k = Tr (T* fc ) M , (B.l) 

where we stress that the trace is over reduced states. Also, let Zj (resp. Z*) be the annulus partition 
function constrained to have exactly 2j unblobbed non contractible loops (resp. 2j — 1 unblobbed and 
1 blobbed non contractible loops). In other words, Zj (resp. Z*) consists of the terms in the full 
partition function Zj<) t M{x,y,l,fri) whose l,m dependence is P J (resp. l 2 ^~ l m). The goal is to search 
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Figure 19: Construction of invariant restricted states, (a) A configuration contributing to Zj with 
N = 12 and j = 2, here depicted as a state, (b) Application on the bottom of the reduced state 
corresponding to the top half of (a), (c) After removal of arcs one has simply 2j links. 



for a decomposition of the form 



Zi = 



z* = 



N/2 



^[DjPVKu + Dgl^mKl 



k=j 

N/2 
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k=j 



D^*l 2k K k +D^l 
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mKl 



(B.2) 



where , D^* , and D^J* are coefficients to be determined. In the notation of (j3.11|) we have 

then 



n* — 
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2.; 
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i=0 
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(i)]2i 



D$*l 2i - l m 



(B.3) 



with D, 



(J ) 



D, 



D 



0> 



0. 



o* 0* " j 

Rather than solving directly for the decomposition of Zj in terms of K k , the idea [26] is now to 
turn the problem upside down and look for the decomposition of K k in terms of Zj : 



N/2-1 N/2 

J t 2j ^ J* pj-l m 

j=k 3=k+l 



N/2-1 



K t = e ^ri+E^ 



J= k 



N/2 



j=k 



3* 



(B.4) 



The determination of the coefficients E can be turned into a combinatorial counting problem as 
follows. First, recall that K k and Kt were defined as traces over restricted states (in contradistinction to 
the partition function which, as we have seen in (|A.2[) . is more complicated than just a trace). We must 
now determine how many times each Zj and Z* occurs within a given trace. Consider therefore some 
configuration C on the annulus that contributes to (say) Zj (i.e., has 2j non contractible unblobbed 
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loops). An example with j = 2 and N = 12 is shown in Fig. 119b . It is convenient to not represent 
the contractible loops within the configuration, i.e., to depict it as a state. This configuration will 
contribute to the trace only over such restricted states S that are left invariant by the action of the 
configuration. Therefore, S must contain the same arcs as does C in its top row (see Fig. [19b). It 
suffices therefore to determine the parts of S which connect onto the starting points of the 2j non 
contractible loops (see Fig. [TOb). Since the goal is to determine the contribution to (say) Kk, precisely 
2k strings and j — k arcs must be used. 

With this in mind, the coefficients E are then determined as the following counting problems. In 
all cases, construct a reduced state on 2j strands, using 2k strings and j — k arcs. Further: 

(k) 

• For E - , all strings are unblobbed, but the exterior arcs to the left of the first string may be 
blobbed. 

(k)* 

• For Ej , the problem is the same, except that the leftmost string must be blobbed. But 

(k)* (k) 

evidently this leads to the same counting, and so E - = E) ' . 

• For Ej^\ the leftmost strand becomes blobbed (since we are considering a contribution from Z*), 
and so must connect onto a blobbed object (arc or string). But as the strings are unblobbed 
(since we are considering a contribution from Kk), it follows that the leftmost object is a blobbed 
arc. 

• For E\ , the leftmost object (arc or string) as well as the leftmost string must be blobbed. 

These counting problems are easily solved using generating function techniques. As a warmup, 
consider the counting of restricted states made up of only arcs. Associate to each pair of sites an 
activity z. A state is either empty, or has a leftmost arc which divides the space into two parts (inside 
the arc and to its right) each of which can accommodate an independent arc state. The generating 
function f(z) therefore satisfies f(z) = 1 + z[f(z)] 2 with regular solution 



m= 1 —^^ = J2c n z-. (B.5) 

71=0 

The coefficients are the celebrated Catalan numbers C n — n \fn+iy. • 

Consider next states made up of only arcs, but in which exterior arcs may (but need not) be 
blobbed. Call the generating function g(z). If the state is non-empty, the leftmost arc is necessarily 
exterior. Inside it are f(z) states, and to its right g(z) states. Thus, g(z) = 1 + 2zf(z)g(z), or 

v n— v / 

We can now attack the case of £f\ Since there are 2k strings (all unblobbed), all of which divide 
the space into independent parts, the generating function reads 

h k (z)=g(z) x z k [f(z)] 2k = E( T1 „ n J z " ( R? ) 

n—k 



and we infer that E ( f } = Ef ] * = (.^ fc ). 

For EfJ it follows fr om the above observations that the generating function with 2k strings reads 



i k (z) = zf(z)h k (z)= L_ k _ 1 ) zn ( R 



n=k+l 



and so = ( .^_\) 



2G 



(k)* 

Finally, for E^J we must distinguish between the cases where the leftmost object is an arc or a 
string. This gives the generating function 



Jk (z) = (1 + zf(z)g{z)) z k [f(z)] 2k = £ ( 2U _ X ) z" 



(B.9) 



whence Ef> = (^) . 

To summarize, we have shown that 
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It is easily checked that this system of equations is invertible, and after some straightforward manipu- 
lations the Zj and Z* can be isolated. In the notation of (|B.2j) the solution reads 



C?) 



(_l)i+fc 



2A- 



2fc 



= (-1 



2k- 1 

, , ( 3 + k - 1 
2k- 1 



(B.ll) 



Applying (1B.3|) then finally leads to (|3Tfl) . 
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